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^SJ ■ rections, one important ingredient is the analytical integration over the one-particle phase 

Tjj- ■ space of the unresolved particle that is necessary when adding the subtraction terms to 

, the virtual corrections. We present the implementation of these integrated dipoles in the 

MadGraph framework. The result is a package that allows an automated calculation for 
the NLO real emission parts of an arbitrary process. 
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1. Introduction 

Multi-particle final states are the basis of many physics studies at the CERN LHC. In 
searching for physics beyond the standard model, one is aiming to identify new particles 
from their decay products, which could often result from decay chains. Likewise, accompa- 
nying final state particles may help to improve the ratio of signal to background processes, 
as done for example in the Higgs boson search through the vector boson fusion channel. 
Meaningful searches for these signatures require not only a very good anticipation of the 
expected signal, but also of all standard model background processes which could result in 
identical final state signatures. From the theoretical point of view, high precision implies 
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that one has to go beyond the leading order in perturbation theory to be able to keep up 
with the precision of the measurements. 

For leading order processes there have been many developments concerning event gen- 
eration and simulation tools in the last two decades such as MadGraph/MadEvent [1-3] 
CompHEP/CalcHEP [4]/ [5], SHERPA [6,7] and WHIZARD [8] and also programs using 
different approaches such as ALPGEN [9] and HELAC [10]. All these programs are multi- 
purpose event generator tools, which are able to compute any process (up to technical 
restrictions in the multiplicity) within the standard model, or within alternative theories 
specified by their interaction Lagrangian or Feynman rules. They usually provide event 
information which can be passed into parton shower, hadronization and/or detector simu- 
lation through standard interfaces. 

Next-to-leading order (NLO) calculations are at present performed on a process-by- 
process basis. The widely-used programs MCFM [11,12], NLOJET+-F [13], MC@NLO 
[14-17] and programs based on the POWHEG method [18-23] collect a variety of differ- 
ent processes in a standardized framework, the latter two methods also match the NLO 
calculation onto a parton shower. The POWHEG box [24] provides a toolkit for adapting 
further NLO calculations to match onto parton showers. 

The NLO QCD corrections to a given process with a n-parton final state receive two 
types of contributions: the one-loop virtual correction to the (2 — t- n)-parton scattering 
process, and the real emission correction from all possible (2 — > n -|- l)-parton scattering 
processes. For the numerical evaluation, one has to be able to compute both types of 
contributions separately. 

The computation of one-loop corrections to multi-particle scattering amplitudes was 
performed on a case-by-case basis up to now, the calculational complexity increased consid- 
erably with increasing number of external partons. Since only a limited number of one-loop 
integrals can appear in the final result [25-27] , the calculation of one-loop corrections can be 
reformulated as the determination of the coefficients of these basis integrals, plus potential 
rational terms. Based on this observation, a variety of methods for the systematic deter- 
mination of the one-loop integral coefficients and the rational terms have been formulated, 
and first fully automated programs for the calculation of one-loop multi-parton amplitudes 
are becoming available with the packages Cut Tools [28,29], BlackHat [30], Rocket [31] and 
GOLEM [32], as well as independent libraries [33]. 

These recent technical advances have allowed the calculation of NLO corrections to 
several 2 — )■ 4 processes, which is the current frontier in complexity. The first result 
in this class of processes were the electroweak corrections to four-fermion production in 
electron-positron annihilation [34], based on [33]. Most recently NLO QCD corrections 
were obtained for genuine 2 — > 4 processes at hadron colliders: H^-|-3 jet production [35-38], 
Z -I- 3 jet production [39], pp — > ttbb [40-42] and pp — > tijj [43], as well as for the quark- 
antiquark contribution to pp — > bbbb [44] . 

The real emission corrections contain soft and coUinear singularities, which become 
explicit only after integration over the appropriate real radiation phase space yielding a 
hard n-parton final state. They are canceled by the singularities from the virtual one- 
loop contributions, thus yielding a finite NLO correction. To systematically extract the 
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real radiation singularities from arbitrary processes, a variety of methods, based either on 
phase-space slicing [45] or on the introduction of process-independent subtraction terms [46] 
have been proposed. Several different algorithms to derive subtraction terms are available: 
residue (or FKS) subtraction [47] and variants thereof [48], dipole subtraction [49,50] and 
antenna subtraction [51-54]. 

Especially the dipole subtraction formalism, which provides local subtraction terms 
for all possible initial and final state configurations [49] and allows to account for radiation 
off massive partons [50], is used very widely in NLO calculations. The generation of 
dipole terms for subtracting the singular behaviour from the real radiation subprocesses has 
been automated in various event generators: in the SHERPA framework [55], the TeVJet 
framework [56], the HELAC framework [57] and in the form of independent libraries [58] 
interfaced to MadGraph. The MadDipole package [59] provides an implementation within 
MadGraph. An implementation of the residue subtraction method is also available within 
MadGraph [60]. 

For a full NLO calculation, the dipole terms have to be integrated over the dipole 
phase space, and added with the virtual corrections to obtain the cancellation of infrared 
singularities. So far, only one of the implementations [57] provides these integrated dipole 
terms including all masses and possible phase-space restrictions, and constructs the corre- 
sponding integrated subtraction terms. It is the purpose of this paper to implement the 
generation of the integrated subtraction terms into MadDipole, which will consequently 
allow to carry out the full dipole subtraction within the MadGraph/MadEvent framework. 
The output results are Fortran subroutines which return the squared amplitude for all 
possible unintegrated and integrated dipole configurations in the usual MadGraph style. 

With a complete treatment of the real NLO radiation, MadDipole is a crucial building 
block of automated NLO calculations. This automation is a very high priority for LHC 
predictions [61], and can likely be accomplished only in a collaborative effort, with different 
groups supplying different pieces, linked through standard interfaces [62]. 

The paper is structured as follows: in Section |2|, we briefiy review the structure of un- 
integrated and integrated dipole terms, Section |3| discusses the expansion of the integrated 
dipoles, and normalization conventions used in this. The MadDipole implementation of 
the integrated dipoles is described in Section |^, with functionality checks documented in 
Section ^. Finally, we conclude with Section ^. 

2. Structure of NLO dipole subtraction terms 

The fundamental building blocks of the subtraction terms in the dipole formalism [49,50] 
are dipole splitting functions ^ij,k, which involve only three partons: emitter i, unresolved 
parton j, spectator k. A dipole splitting function accounts for the collinear limit of j with 
i, and for part of the soft limit of j in between i and k. The dipole factors, which constitute 
the subtraction terms, are obtained by multiplication with reduced matrix elements, where 
partons i, j and k are replaced by recombined pseudo-partons ij, k. The full soft behavior 
is recovered after summing all dipole factors. 
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Throughout the whole paper we are using the notation introduced in Refs. [49] and [50]. 
The dipole factors are subtracted from the real radiation contribution at NLO. They sub- 
tract the singular contributions where one parton from the NLO real radiation contribution 
becomes soft and/or collinear, such that the phase-space integral of this contribution can be 
preformed numerically, including arbitrary kinematic restrictions on the final-state phase 
space. 

The dipole factors are integrated analytically over the dipole phase space (which fully 
includes the infrared singular regions), such that the integrated dipole factors have the 
same kinematic structure as the virtual one-loop NLO corrections and the collinear coun- 
tertcrms from mass factorization. These integrated dipole contributions can then be added 
to the virtual and mass-factorization corrections, thereby accomplishing the cancellation 
of infrared singularities. 

Several algorithms to automatically generate the unintegrated dipole terms for arbi- 
trary processes have been devised and implemented in various matrix element generator 
frameworks. They are available for SHERPA [55], HELAC [57], for MadGraph [59] as well 
as a implementation of stand-alone routines [56,58]. Also implementations of integrated 
dipoles are available in the same codes. However, only the implementation based on the 
HELAC framework has the full set of integrated dipoles for arbitrarly masses and phase- 
space restrictions [57]. It was used in the context of the calculations of NLO corrections 
to pp — > ttbb [42] and pp ttjj [43]. In this work, we document our implementation 
of the integrated dipoles in MadDipole, which is a package in the MadGraph framework. 
With this extension, MadDipole computes the full NLO dipole subtraction of real radia- 
tion and performs the infrared cancellations in a fully automated manner for color- and 
helicity-summed matrix elements squared. Our implementation was already used in the 
computation of the NLO corrections to pp bbbb in the quark-initiated channel [44] . 



2.1 Phase space restriction: a-parameter 



The calculation of the subtraction terms is only necessary in the vicinity of a soft and/or 
collinear limit. Away from these limits the amplitude is finite and there is in principle no 
need to calculate the computationally heavy subtraction terms. The distinction between 
regions near to a singularity and regions without need for a subtraction can be parametrized 
by a parameter usually labelled a with a G [0,1], which was introduced in Ref. [63] for 
processes involving partons only in the final state. The case with incoming hadrons, i.e., 
with partons in the initial state, is described in Ref. [13]. 

Using the notation of Ref. [13], the contribution from the subtraction term to the 
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differential cross section in the real radiation channel can be written as 



{n+1} '^{"+1> 



:< XI 'I^ijAPa,Pb,Pl,---,Pn+l)F!f\pa,Pb,Pl,--,Pij,Pk,--,Pn+l)&{yij,k < 

{, pairs ky^ij 



pairs 



T^ij{Pa,Pb,Pl, ■■■,Pn+l)Fj'\pa,Pb,Pl, ■■,Pij, ■•,Pn+l)0(l - Xij^a < a) 



+{a -H- b) 

+ [^kiPa,Pb,Pl, ■■■,Pn+l)Fj'^ ipa,Pb,Pl, ■■,Pk, ■■,Pn+l)Q{Ui < a) + (a O 5) 
+ Y1 [T^''''\Pa,Pb,Pl,---,Pn+l)F!f\pa,Pb,PU---,Pn+l)&iVi < a) + (o O 6)] | . 

(2.1) 

The functions Pjj.fc, T^fj, P^* and D"*'^ are the dipole terms for the various combinations 
for emitter and spectator. ^^^^+1} denotes the summation over all possible configurations 
for this (n + l)-particle phase space which is labelled as (ir("+^) and the factor S^n_^_iy is the 
symmetry factor for identical particles. In MadDipole, we have introduced four different 
a-parameters, one for each type of dipoles [59]. In our code they are called alpha Jf, 
alpha_f i, alpha_if and alpha_ii for the final-final, finial-initial, initial-final and initial- 
initial dipoles, respectively. The actual values for these parameters are by default set to 
unity, corresponding to the original formulation of the dipole subtraction method [49,50], 
but can be changed by the user. 

The integrated dipole factors, which are to be added with the virtual n-parton contri- 
bution, also depend on a. For case of massless partons, the a-dependence of the integrated 
terms is stated in [13, 63] while for massive partons results for most cases can be found 
in [57,64,65]. The remaining cases, i.e., the (finite) massless-to-massive splittings, can be 
found in the appendix. 



2.2 Regularization scheme dependence 

Calculating objects that contain divergences requires a systematic prescription of how to 
deal with these divergences, i.e., a scheme for their regularization. In NLO calculations, the 
same regularization scheme has to be applied in the real emission part and in the virtual 
corrections. Both these contributions will differ between different regularization schemes, 
while their sum (i.e., the full NLO result) is scheme-independent. Therefore it is necessary 
to clearly specify which regularization scheme one is using. 

In QCD calculations, there arc mainly two types of regularization schemes used, namely 
dimensional regularization [66-69] and dimensional reduction [70-72]. Both extend the 
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dimensionality of space-time to d = 4 — 2e, resulting in divergences becoming explicit as 
poles in 1/e. A discussion about their subtypes and their differences can be found in [73]. 

In the real radiation contribution, the dependence on the regularization scheme does 
not yet appear explicitly at the level of the unintegrated dipole terms, and we consequently 
did not address this issue in the previous release of MadDipole [59] . 

The helicity subroutines on which MadGraph and MadDipole are build evaluate matrix 
elements in four dimensions. Therefore we can compute the subtraction terms only in 
regularization schemes in which the external particles are defined in four dimensions. The 
two most widely used are the 't Hooft-Veltman scheme (tHV) in dimensional regularization, 
which is the default used in our implementation, and the four-dimensional helicity scheme 
(FDH). Both methods differ only by a finite shift [49,74]: 

V,(e)*HV ^ ^FDH(^) ^ ^^(^)tHV ^ (2.2) 

lq = lq = \cF , % = ^Ca . (2.3) 

In the massive case there is no dependence on the regularization scheme [75]. There is a 
simple flag in our code that allows to change between these two schemes. 



3. Expansion of integrated dipoles in e 



Integration of the dipoles makes their infrared singularities explicit as poles in the di- 
mensional regularization parameter e. The formal structure of the integrated dipoles is 
independent of the configuration we are considering. For definiteness, we discuss only the 
final-final case, the same structure also holds in all other cases. For initial state hadrons, a 
additional collinear contribution is present, which is rendered finite by mass factorization, 
which we will describe here as well. 



3.1 Infrared final-state singularities 

In the final-final case, the integrated dipole function is written as 



/ 



{Pi+Pj)-mf- 27rT{l - €)\2pij -pk,' 



Depending on the configuration, the splitting function and the propagator on the left hand 
side of (3T) change their form; the structure of the result on the right hand side remains 
the same. 

The factor Vij{e) is determined by the specific configuration. It is singular in the limit 
e — > 0, and is expanded as a Laurent series in e. Prefactors taken out from the expansion 
must be consistent between the dipoles and the virtual one-loop corrections to the process 
under consideration. 
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In the MS-scheme, only universal factors are taken out, and this expansion can be 
written symbolically as 



where 7 is the Euler-Mascheroni constant, 7 = 0.5772 .... This structure is the basis for 
our implementation. 

The specific values of Uij-i, yij-2 and yij-s depend on the splitting one is considering. For 
instance if one takes the massless final-final quark-gluon splitting, i.e., where the emitter 
is a massless quark int the final state, the unresolved particle is a gluon, and the spectator 
is also a massless final state particle, then the integrated splitting function is given by 



1_ 3_ _ 7r2 
e"2 + 2^ 



Vqg{e) = CF(- + - + 5- — ]. (3.3) 



Consequently, the expansion coefficients in ( |3.2D become 



27r V2 \2pij-pk 

as „ / 1 , 2 /i^ \ 3 - / fi^ \ Itt"^ 



W = ^CF^log^ +^log -1- . (3.4) 

27r \2 \2pij-pkJ 2 \2pij-pkJ 12 J 

3.2 Initial-state collinear behaviour 

The cases with initial state radiation, i.e., initial-initial and initial-final, are slightly more 
involved. Not all singularities that occur in the real emission process are cancelled by the 
virtual corrections. The integrated initial-final dipole functions take the form: 



/ 



[Pi +Pj) - m-j ns{a) 
- 2vrr(l- 6) Upa-pJ ^ ' ' 

= [-^ + + ) ■ (3.5) 

The left-over singularities arise from collinear splitting off the initial emitter particle. For 
example, for the initial-final dipole describing the gluon radiation off an incoming quark, 
one has 



1 



V<?'9(z;e) = —Pii{z)+6{l-z) 



€ 



2 



V,,(e) + (^-5)CF 



+ 55'9(^) + 0(e), (3.6) 



where B'^''^{z) contains regular functions and plus-distributions in z. 
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The collinear singularity remains and is absorbed into the parton distribution function. 
This is done by introducing a coUinear counterterm and its contribution to the cross section 
is given by (6.6) of [49]: 

da^ip;^,^) = ^^J—.^. I' dz [-1 (^^y p-\z) + K^'liz)'^ da^{zp) , (3.7) 

where in MS scheme, Kp'^{z) = 0. 

Neglecting the sum over the partonic subprocesses, we have a counterterm contribution 
of the form 



27rr(l-e) V 

which is added to the integrated dipole. 

In the MS scheme the expansion of (|3.8|) is given by 



2\ e 



(3.8) 



The coefficients and ^.3 are given by 



<« = £-^~'(=)log(|j^ (3^10) 
3.3 Structure of program output 

The final output the program provided to the user is then given by all contributions of 
the coefficients yi and yf respectively multiplied with a born level matrix element that is 
modified by its color structure. This explicitly means 

js- Vi ■ m{'^---m\ L,i '' \l---m)rn , 

\- (^2+^2)- , 

T T 

finite: {y-i + vl) ■ .m(l • • • rnl— ■ ■ ■ m)^ . 

i 

The contributions from the collinear counterterms are of course only present if there are 
initial state QCD particles involved. By adding the collinear counterterm, only endpoint 
singularities, which occur at z = 1, remain. Those then completely cancel with the virtual 
corrections. The finite pieces can contain regular functions of z as well as 5-functions and 
plus-distributions. 

4. Implementation and how to use the integrated dipoles 

The installation and running of the new package is very similar to the already existing 
MadDipole package: 
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1. Download the MadDipolc package (version 4.4.35 or later), MGjyiE_DIP_V4 . 4 . ?? . tar .gz, 
from one of the MadGraph websites, e.g., http://madgraph.phys.ucl.ac.be/. 

2. Extract and run make in the MadGraphll directory. 

3. Copy the Template directory into a new directory, e.g., MyProcDir to ensure that 
you always have a clean copy of the Template directory. 

4. Go to the new MyProcDir directory and specify your process in the file . /Cards /pro c_card . dat. 
This is the (n + l)-particle process you require the subtraction term for. 

5. Running . /bin/newprocess generates the code for the (n + l)-particle matrix el- 
ement and for all dipole terms and their integrated versions. After running this 
you will find a newly generated directory . /SubProcesses/PO_yourprocess {e.g., 
./SubProcesses/PO_e+e-_uuxg) which contains all required files. 

In the . /SubProcesses/PO_yourprocess directory all the files relevant to that par- 
ticular subprocesses are generated. In particular this includes the {n + 1) particle matrix 
elements in the file matrix. f and the dipoles in the files dipol???.f, where ??? stands 
for a number starting from 001. Furthermore the directory has two files, dipolsum.f and 
intdipoles . f , where the sum of the dipoles and their integrated versions are calculated, 
respectively. 

Here we discuss in more detail the syntax and implementation of the integrated dipoles 
in the intdipoles . f file. In this file there are two subroutines of the form 

intdipoles (P, X, Z, PSWGT EPSSQ, EPS, FIN) and 
intdipolesfiniteCP, X, Z, PSWGT EPSSQ, EPS, FIN), 

where the input parameters are the phase space point P(0:3,nexternal), the Bjorken 
X values of both incoming parton distributions, X(l:2), and the momentum fraction of 
the incoming parton that goes into the hard process after an initial state radiation or if 
the spectator is an initial state particle, Z. It is the latter quantity which is denoted with 
Xij,a, Xik a a-nd Xi^ab respectively in [49], but wc shall refer to it here simply as z. In the 
numerical integration, it must be taken uniformly over the interval z € [0; 1] to ensure 
correct representation of the distributions. Furthermore, the phase space weight should 
be passed as well. These first four arguments should be provided by the user. For the 
given phase space point, the routines evaluate the sum of all integrated dipole subtraction 
terms (integrated dipole factors multiplied with the appropriate reduced matrix elements) 
after mass factorization of the collinear singularities. The routines call external subroutines 
(explained below) from dipolesum. f which supply the parton distributions appearing with 
the reduced matrix elements, apply event rejection cuts and pass the event information into 
histograms. 

The integrated dipoles in the two subroutines correspond to the unintegrated dipoles 
in dipolsumC . . ) and dipolsumf inite ( . . ) . The dipoles in the dipolsumf inite ( . . ) and 
intdipolesf inite ( . .) subroutines are not needed to cancel singularities because they 
correspond to gluon splittings into massive particles. However, they can be useful for 
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checks when taking the hmit of vanishing quark masses, or to cancel some large logarithms 
in the real emission matrix elements [59]. 

The output parameters are the coefficients of the divergent and finite terms: EPSSQ 
is the coefficient of EPS is the coefficient of 1/e. After inclusion of the collinear 

counterterms they contribute with a factor 6{1 — z), and are real numbers. FIN is the 
finite coefficient. The calculation of the FIN coefficient requires some explanation, since 
its coming from distributions in z. In the intdipoles ( . .) and intdipolesf inite ( . .) 
subroutines a three-dimensional array FINITE is filled with the contributions from the 
various dipoles: 



FINITE(l) 
FINITE(2) 
FINITEO) 



regular function in z, 
coefficient of 5(1 — z), 
coefficient of 5(z+ — z). 



The last entry appears only for massive dipoles, with 

z+ = l 



(4.1) 



where hq is the rescaled fermion mass occurring in the splitting. 

In this decomposition, the (5-functions and (+)-distributions are carried out by taking 
into account the convolution with the product of a reduced matrix element g{z) (generated 
by MadDipole out of MadGraph) and a parton distribution function h{z) (supplied by the 
user through the subroutine dipolepdf ) 



dz6{l- z)g{z)h{z) 



dz{f{z)),g{z)h{z) 



dz g{l)h{l) 



in FINITE(2) 



dz 



f{z)g{z)h{z)-f{z)g{l)h{l) 



(4.2) 
(4.3) 



in FINITE(l) 



in FINITE (2) 



In the massive case the point z = 1 can not be reached in all cases but we may have a 
reduced endpoint z+, such that instead of a 5(1 — z) we then have a 5-distribution of the 
form 5(z+ — z) which is the third entry of the array of the finite terms. If we have such 
a reduced endpoint then also the (+)-distribution is generalized into a 2:+-distribution as 



defined in ( A. 13 ). As before, we have the following implementation: 



dz5{\ — zj^)g{z)h{z) 



dz g{z+)h{z^ 



(4.4) 



in FINITE (3) 



/" 

Jo 



dz{f{z))^^g{z)h{z 



dz@{z+-z) f{z)g{z)h{z)-f{z)g{z+)h{z+) . (4.5) 



in FINITE(l) 



in FINITEO) 



By making a transformation of variables and computing the parton distribution func- 
tion not at X (Bjorken's x) but at x/z, the matrix element itself becomes independent 
of z as the initial state particle, that radiates the unresolved particle, then comes with 
the momentum fraction of x/z ■ z = x. Therefore the set of momenta which are used to 
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calculate the matrix element do not depend on z. This means that g{z) 
Eqs. (pl-i^D. 

The final result FIN is the sum of the three contributions: 



= 5'(^+) = 5(1) in 



FIN = (FINITE(I) jz + FINITE(2) + FINITE(3) /z+) x PSWGT. (4.6) 

Since FINITE (3) appears only in massless-to-massive splittings, it is always zero for the 
intdipolesC . . ) subroutine. Conversely, the EPSSQ and EPS should be zero when coming 
from the subroutine intdipolesf inite ( . . ). 

For ease of use, we have provided three dummy subroutines/functions that the user 
might want to fill when using the code. In the code there are consistent calls to these 
subroutines. These subroutines/functions can be found in the file dipolsum. f and are 

passcutsdip(P): In this LOGICAL FUNCTION the user should provide a set of cuts that he 
wants to be applied to the phase space points. Note that in general the phase-space 
mapping for each unintegrated dipole is different; there needs to be a call to this 
function for each dipole. It should return FALSE if the point fails the cuts. By default 
every point passes the cuts. 

dipolepdf (P,legl,leg2,WGT): In this SUBROUTINE the user should provide the value for 
his/her favourite PDF set for the two incoming particles with PDG codes passed 
by legl and leg2. The factorization scale should be defined in the include file 
dipole. inc. The weight from the PDF should be returned in the argument WGT, 
which is set to 1 by default. 

writehist(P,WGT): In this SUBROUTINE the phase-space point (for each dipole) are pro- 
vided together with its weight. The user could use these to fill histograms or save 
ntuples. 

Besides these three subroutines/functions, the more technical parameters, like the q- 
parameter (see Sec. |2.l| ), the number of fiavors, the renormalization scheme and the scales 
can be set in the include file dipole . inc. When changing any of the parameters in this file 
the code should be recompiled (after removing the object files) for these changes to take 
effect. 

Besides the already existing check checking program, that checks the limits of the real 
emission matrix element minus the subtraction terms, we provide the user with another 
sample program, checkint, to calculate the value of the integrated subtraction terms for 
a given (or random) phase space point. 

More details and latest news, updates, bug fixes, etc. can be found at 

http : / / cp3wks05 . f ynu.ucl . ac .be/twiki/bin/view/Sof tware/MadDipole. 
5. Checks 

To verify the implementation we have performed two different kinds of checks: indepen- 
dence on the phase space restriction parameter a and comparison with the implementation 
of dipoles in the MCFM program [11,12]. 
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(a) (b) 
Figure 1: Two examples for the cancellation of the a-dependence between unintegrated and 
integrated dipoles on specific contributions. Since these do not correspond to full physical processes, 
we use an arbitrary normalization. The plot on the left shows the a dependence of the initial-initial 
dipole where a gluon is radiated of the initial quark lines. On the right the a dependence of an 
initial-final dipole where the spectator is massive is shown. For both figures only the a parameter 



shown is varied, all others are kept fixed. The upper dashed line is the first contribution of (5.1), 
i.e., the sum of matrix element and unintegrated dipoles. The lower dashed line is the second 
contribution, i.e., the finite terms of the integrated dipoles. The (red) solid line is the sum of both. 
For the sum we also included the Monte- Carlo error to show that the results for different values are 
consistent with each other. 



Both the non-integrated and the integrated dipoles depend on the a-parameter, the 
dependence on this parameter should cancel in the sum. To validate this, we require: 

/ (d(T^ — da^) + / (finite parts of int. dip.) = const, (5.1) 

which must be a constant in that sense that it should not depend on a. We have vahdated 
this for all 27 different cases (emitter /spectator in initial/final state, and mass assign- 
ments) listed in the appendix. Figure |l] shows the dependence on the a-parameter for two 
examples: in the left plot is the dependence shown on the ajj-parameter that governs the 
initial-initial dipole phase-space restriction, and in the right plot for a processes with mas- 
sive spectators the dependence on the a, /-parameter (that restricts the initial-final dipole 
phase space) is shown. The solid (red) lines display the quantity defined in (^), which is 
observed to be independent on a over several orders of magnitude. 

This is a very powerful check because it includes several aspects of the implementation. 
It verifies not only the correct implementation of the ^-functions for the unintegrated 
subtraction terms and the a-dependent correction terms for the integrated dipoles, but 
independence on this parameter can only be achieved if the correct parton distribution 
functions are called with the right arguments, and if the various terms contributing to 
the integrated dipoles are summed correctly. Moreover, an inconsistent infrared-unsafe 
implementation of the cuts (by the user) will lead to a dependence on the a parameter. 

The a-independence does, however, not probe important features related to the pole 
structure of the subtraction terms, namely the scheme dependence and the factorization 
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scale dependence, neither does it show other finite contributions. Therefore our second 
check was a direct comparison of our implementation against MCFM [11, 12]. We have 
compared the results for single phase space points, which allowed us to directly probe 
the implementation of the integrated terms. The various terms (regular functions, S- 
and plus-distributions) could be checked separately. We note that not all possibilities 
(massless/massive, initial/final emitter/spectator) are present in MCFM and therefore not 
all could be checked. This includes in particular the finite dipoles: massless-to-massive 
splittings have no divergences, but do have a potentionally large logarithm which might be 
useful to subtract. 

6. Conclusions 

The MadDipole package [59] provides dipole subtraction terms for the evaluation of real ra- 
diation corrections in NLO calculations within the framework of the MadGraph/MadEvent 
matrix element and event generator [1-3]. In this work, we described the extension of the 
MadDipole package to include the integrated dipole terms, which are required to combine 
the real radiation corrections to a given process with the virtual corrections and coUinear 
counterterms of the parton distributions. With the newly developed subroutines, Mad- 
Dipole provides the unintegrated and integrated dipole subtraction terms. 

The integrated subtraction terms are convoluted with the user-supplied parton distri- 
butions and are combined with the coUinear counterterms from mass factorization. Con- 
sequently, infrared singularities appear only in the kinematic endpoints, which correspond 
to the kinematics of the virtual corrections. The MadDipole output can thus be readily 
combined with the results from one-loop matrix element generators, which are currently 
under rapid development [28-32]. 

A first application of the MadDipole package, in combination with the GOLEM one- 
loop amplitude generator [32], was in the calculation of NLO corrections to the quark- 
antiquark contribution to pp bbbb [44] . Many more applications are likely to follow. 
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A. a-dependence of the integrated dipoles 

In this appendix we give a list of the references, where the ct-dcpcndcnt terms can be found 
in the literature. For the few cases which were not known previously, we give the result 
below. 
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A.l Final-final 



Configuration 



Reference 



Remarks 



(a) 

(b) 

(c) 

(d) 

(e) 

(f) 

(g) 
(h) 

(i) 
(j) 



Q 
Q 
q 
q 

9 
9 

9 - 
9 - 

9 
9 



> Q g, nik > 
>Q g, mk = 

> q g, mfc > 

> q g, rrifc = 

> 9 9, rnk> 

> 9 9, "ifc = 
Q Q, ruk > 
Q Q, ruk = 

> q q, ink > 

> q q, mfc = 



[42] (A.9,A.17,A.20) 
[64] (A.22,A.23) 
[64] (A.29,A.31) 
[63] (22) 

[42] (A.9,A.19,A.22) 

[63] (24) 



Different definition of a 



this work, ( A. 6 ) 
this work, (JjA^ 
[42] (A.18,A.21) 
[63] (23) 



Different definition of a 



Different definition of a 



In the massive case there is an ambiguity how one defines a. The upper limit of the 
integration over y is given by the variable y^. One can now define a in such a way that 
the maximal value for a is given by y+. In that case the integration range is given by 
the interval y € [«,?/+]. Another possibility is to have a = 1 as the maximal value which 
implies that the allowed integration range is given by 

ay+ < yij,k < y+- (A.i) 

The first definition has been used in [42] whereas we use the latter. 

The two cases (g) and (h) are in principle not needed as the result is still finite due 
to the masses of the quarks. For the sake of completeness we also included this cases and 
therefore give the result. 

A.1.1 Case (g) 

Case (g) is the splitting of a gluon into two massive quarks {g ^ Q Q) with a massive 
spectator. The splitting function for that process is given in Eq.(5.18) of [50], 



(VQQ,fc) = 8V^asrR— <;i 



Zi{l - Zi) - (1 - k)z+z_ 

(A.2) 



2^^ + (1 - 2/.^ - fii)y,,^ 

where we neglect the last term because of k = in our implementation. The integration 
ranges are 

1 ~ Vij iVij h , 1 + Vij iVij h , . , 

— — < Zi < — —. (A. 4) 
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Here, the integration over j/jj ^ does not start at /j = but at a value larger than 
zero. But introducing the a parameter and calculating the correction terms implies new 
integration boundaries for the yij^k integration according to ( A.l): 



^ -2 < ay+, (A.5) 



which implies that a must not be chosen to be too small for this splitting. With this 
restriction, we find 



^li °) = (^) - Tk ( (bd (-8aM,' log (^ac'y+ - - 4^i* - 4^^) + 2ac^ log (^ac^y+ ~ d^ - 4/i| - 4^*) 



+2aclog ya.(?y+ — — 4/i| — 4/i^j + 4a^^ log ^ac^y+ — dy — 4/i| — 4/i* 

-2a (cV c - 4^t| + 2^1) log (^(1 - Qy+) (-c - 2/^1)'^' (c - 2^| + l)) 
+8a/i;* log {-b^c^ - 4^*4 + c''y+ - 4^i|) - 2ac'' log {-h^ - 4^^ + c'y+ - 4^t*) 



-2aclog (-b^c2 - 4^ij + c^y+ - 4^i|j - ^afi] log (-b^c^ - 4^** + c^y+ - 4^** 
+2a (c^ + c - 4^* + 2^,2) log ({y+ - 1) (-c - 2^*,^)'/') (c - 2^,] + l)) 
-30^2^1 - clog(-2(acy+ + d)) + 4cm- y'2/i2 - clog(-2(acy+ + d)) 
-2cy^2/i2 -clog(-2(acy+ + d)) + 30^2^*2 _ clog(-2(6 + cy+)) 
-4c/i,y2M|-clog(-2(& + cy+)) + 2c^2^if _ clog(-2(6 + cy+))) 

2/i| \ 



^-Ihd^lp.) - c {c - 2(c + + 4^i*) tan~^ 



a^c^y\ - 4^1 



+c 



^2m?-c (c2y+ - 2a6 - d(y+ - 2)) + Ac^,]{h{ay+ - 1) + d(-y+) + d) + 4^(6 - d)) 



2bdJ2n'j ~c(c - 2(c + + 4^i*) tan"^ 



2m? 



/ (3c (2m| - C) ^C2y2 _4^4 ^„2,2j^2 _ 4^4^ j (A.6) 

where we used the following abbreviations: 



c = -1 + 2^1 + ^1 



A. 1.2 Case (h) 



Case (h) describes also the splitting of a gluon into a massive quark pair {g ^ Q Q) however 
with a massless spectator. While the splitting function is the same ([A.2[ ), the integration 
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boundaries are different, namely 



2 < yij,k < 1 (A.7) 

i 



1 - M 



<z, < ' (Ai 



2 2 

The phase space integral gets multiplied by Q{a — i/ij^k) and integration leads to 



2m? 



-log -2 A/a2 (1 - 2^*2)^ _ 4/i4 + Q (2^*,^ - l) + 2tan-^ 



2/i? 



+ log (-2 (2m| + ^1 - - l) ) - 2tan-^ ( ^== ) ) + a/1 - 4^? | 1 • (A.9) 



A. 2 Final-initial 



Config 


'uration 


Reference 


Remarks 




(a) Q 


g 


[64] (A.13,A.14) 






(b) q 


Q g 


[64](A.17) / [13] (11-16) 


Different approaches 




(c) 5 


^QQ 


this work, ( A..12, A.-141 A..iq) 






(d) 5 


q q 


[13] (11-16) 


Different approach in 


MadDipole 


(e) 9 


-^9 9 


[13] (11-16) 


Different approach in 


MadDipole 



The cases (b), (d), and (e) can be found in [13]. Their result however contains already 
the sum of different contributions making use of the /- and X-flavor kernels. As it mixes 
different contributions this is not suitable for the MadDipole implementation. Therefore 
we use [64] for (b) and derive results for (d) and (e) following the approach in [64]. Of 
course, both approaches are equivalent. Again, for the sake of completeness we also add 
the finite case (c). 

A. 2.1 Case (c) 

The one particle phase space for final-initial dipoles is given in Eq.(5.48) of [50]: 

j [dPiiPij;Pa,x)] = ^(27r)"^+^'(2^ijPa)^~' J dXij^aS{x - Xij^a) (l - X + flfj)'^ 

X / d'^-'^Q [ ^ dzi [z+{x) - z^]-' [zi - z_{x)]-' , (A.IO) 

J Jz.(x) 

and the integrated splitting function is given by Eq.(5.53) of [50] as: 
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For case (c), the splitting of a gluon into massive quarks {g ^ Q Q), the integrated spitting 
function is given in Eq.(5.57) of [50]: 

6) = Tr { [J-q{x, mq)].+ + 6{x+ - X) e) + JJ^'Cmq)] } + 0(e) , (A.12) 

where the X4--distribution is defined as: 



(/(.)) .(.).r../,.,e(., 



(A.13) 



Imposing the cut on the a-parameter imphes that the phase space in ( |A.10| ) is multiplied 
with 0(a — 1 + Xija)- 

This leads to a modification of the x+-distribution terms and we get in analogy to 
Eq.(5.62) of [50] 



2 / 1 - X + 2^2 



Q 



[l-xf 



1 — X 



(A.14) 



where we define: 



dx f{x) (g(x))^_ 



dx g{x) (fix) - /(I)) 



(A.15) 



l-a 



The non-singular terms J°^^^{^q) receive: 



Jqq if^Q^oi) = j;-^ (mq) + - I I -4/iQ 



9 



+ 4 -5 



-16^ + 5) / (^1 - 4^2^) + 61og (^Q - 4^^ + 7^) - 6 log (^1 - 4/.^ + l)) . (A.16) 
A.2.2 Case (d) 

Case (d) is just the limit /ng — )■ of ( A.14| ), i.e., 



1 



l-a 



3 VI -X 

which leads to the following additional non-singular terms: 

ja-.m 



Jqq (0,a) = Jgg (0) + -loga. 



(A.17) 



(A.18) 



A.2.3 Case (e) 

In the case of the splitting {g ^ g g) the general structure of the integrated splitting 
function is given by Eq.(5.66) of [50]: 



I'^ggix; e) = 2Cp^{[J''gg{x)]+ + (^(l " x) ^^^{6)] + 0(e) . 



(A.19) 
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The first term [J^g(x)]+ contains all +-distributions but is not a +-distribution itself. In 
the presence of the a-parameter we find 

[J-(x,a)U=(-^ln-^-^-^) +^ln(2-x)G(a-l + x), (A.20) 

which leads to a modification of the terms proportional to 6{1 — x) of Eq.(5.68) of [50] of 
the following form: 

Jgf (x, a) = J^fix) - log2 a - ^ log a. (A.21) 

A. 3 Initial-final 



Configuration Reference Remarks 

(a) ij : q, emitter : g [64] (A. 9, A. 11)/ [13] (11-16) [13] only massless 

(b) ij : g, emitter : q [65] (A. 8) 

(c) ij : q, emitter : g [65] (A. 10) 

(d) ij : g, emitter : g [65] (A. 11) 



From the analytical point of view the limit of a vanishing spectator mass can be performed 
without any problems. However taking the result for a massive spectator and setting the 
mass to zero in the numerical implementation causes problems for the cases (b) and (d). 
For these two cases we calculated the limit analytically and implemented a massive and a 
massless version. 



A. 4 Initial-initial 





Confi^ 


juration 




Reference Remarks 


(a) 


ij : 


q, 


emitter : 


Q 


[64] (A.4)/ [13] (11-16) 


(b) 


ij : 


9, 


emitter : 


q 


[64] (A.5)/ [13] (11-16) 


(c) 


ij : 


q, 


emitter : 


9 


[65] (A.4)/ [13] (11-16) 


(d) 


ij : 


9, 


emitter : 


a 


[65] (A.5)/ [13] (11-16) 



For our implementation the results from [64,65] are used. 
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